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Abstract 

Instanton effects and three-loop perturbative contributions are incorporated into 
QCD sum-rule analyses of pseudoscalar {J^^ = '') gluonium. Gaussian sum- 
rules are shown to be superior to Laplace sum-rules in optimized predictions for 
pseudoscalar gluonium states in the presence of instanton contributions. The Gaus- 
sian sum-rule analysis yields a pseudoscalar mass of (2.65 it 0.33) GeV and width 
bounded by F < 530 MeV. The Laplace sum-rules provide corroborating evidence 
in support of the 2.7 GeV mass scale. 

Key words: pseudoscalar gluonium, glueballs, QCD sum-rules, instantons, PACS: 
14.40. Cs, ll.55.Hx, 11.40.-q 



1 Introduction 



The gluon self-interaction in QCD suggests the existence of bound gluonic 
states known as gluonium or glueballs. The properties of gluonium states have 
been studied within a wide variety of theoretical methods including lattice 
simulations [1,2], QCD sum-rules [3,4,5,6,7,8,9,10,11], and other phenomeno- 
logical models [12]. From the experimental viewpoint, a number of scalar and 
pseudoscalar isoscalar states exist which are potential gluonium candidates 
[13], but a consensus on the nature of these states has not been achieved. 
Similarly, theoretical investigations of the gluonium mass spectrum continue 
to be refined. 

The initial QCD sum-rule analyses of scalar gluonium [4] and pseudoscalar 
gluonium [5] have been extended to include higher-loop contributions of QCD 
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condensates and higher-loop perturbative effects [6,7]. The non-leading con- 
densate corrections arising from the ghion condensate (aG"^) [14] are particu- 
larly significant since they provide the leading contribution in many sum-rules. 

It is also known that (direct) instantons contribute to the pseudoscalar and 
scalar channels [15]. These contributions have a significant effect on the scalar 
gluonium Laplace [3,8,9,11] and Gaussian sum-rules [16]. In particular, the 
instanton effects improve the self-consistency of the lowest-weighted Laplace 
sum-rule containing the low-energy theorem and the higher-weight sum-rules 
which are independent of the low-energy theorem [9,16]. However, the role of 
instantons in the Laplace sum-rules for pseudoscalar gluonlTim have not been 
studied in similar detail. The self-dual properties of the instanton implies that 
the leading instanton effects in the scalar and pseudoscalar gluonic channels 
differ only by an overall sign. This presents the interesting possibility that 
instantons are responsible for the large mass splitting between scalar and 
pseudoscalar gluonium states observed in lattice simulations [1]. 

In this paper, instanton and higher-loop perturbative effects are incorporated 
into the Laplace and Gaussian QGD sum-rules for pseudoscalar gluonium. In 
Section 3 it is demonstrated that the Laplace sum-rules fail to yield an opti- 
mized mass prediction for the pseudoscalar gluonium state. This shortcoming 
is shown to be overcome by the Gaussian sum-rule approach presented in 
Section 4. The full phenomenological analysis of the Gaussian sum-rules for 
pseudoscalar gluonium is presented in Section 5. 



2 QCD Laplace Sum-Rules for Pseudoscalair Gluonium 

Pseudoscalar gluonium is studied through the two-point correlation func- 
tion of the renormalization-group invariant currents j{x) — asC^^GI^^ — 

n (g') = i I d^a;e'''-^(0|r[J(a;)J(0)]|0), ^ -q^ > . (1) 
The correlation function (1) satisfies the dispersion relation 

oo , , 

n (g^) = n(o) + g^n'(o) + ^g^n"(o) - g^i / dt ^f^^ (2) 

to 

relating the QGD prediction 11 (g^) to the hadronic spectral function p{t) 
appropriate to pseudoscalar gluonic states above the physical threshold . 
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The Laplace sum-rules resulting from (2) are 

1 °° 

A(t) = - / t^e-'^pit) dt , A; = 0, 1, 2 . . . , 



(3) 



where the theoretically-determined quantity jCkir) is obtained by applying the 
Borel-transform operator B 



lim t^lA: 
n q'I^ V{N) \dQ\ 



N 



(4) 



iV/Q2=r 

to the appropriately-weighted correlation function 

1 



A(t) = -B 

T 



(5) 



In the resonance (s) plus continuum model [17], hadronic physics is locally dual 
to the QCD prediction for energies above the continuum threshold t = Sq 



p{t) = p'^^'^t) + e{t- So) Imn'5^^(i) 
The QCD continuum contribution 



CO 

(r, So) = - / t^e-*^Imn«^^(i) dt 

71 J 



(6) 



(7) 



so 



is thus combined with the quantity jCk{T) because both are QCD predictions, 
resulting in the following Laplace sum-rules relating QCD to hadronic physics 
phenomenology: 



so 



TT J 

to 

Rk (t. So) = Ck (r) - Ck (r, so) ■ 



(8) 
(9) 



The field-theoretical (QCD) calculation of 11 (Q^) contains perturbative (loga- 
rithmic) corrections known up to three- loop order in the chiral limit of n/ = 3 
massless quarks in the MS scheme [18], QCD condensate contributions [7,20], 
and direct instantons [3,21] 



n(Q^) =0'iog(^ 



ao + ai log ( ^ 1 02 log^ ( % 



i/2 



+ 



bo + bi log I ^ 



c^G') + co-^ {gG') + do-^, {Os) (10) 



- 327r2g4 J p4 



dn{p) 
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where 



1 + 20.75- + 305.95 - 
TT Vtt 



ai 



\3 rg 



a 



, + 72.531 

4 TT, 
2 



, 02 



10.1250 (^^^ 

2 



6o = 47r- , 6i = 97r - , Cq = -Stt - , do ^ Sir 
TT \7r/ Vtt/ 



TT 



9 f abcG ii^yG pG 



(08) = ( {o^fabcGlpl,) ) + 10 ( (a/a6cG;;,Gj,) 
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(11) 

(12) 

(13) 
(14) 
(15) 



and K2{x) is the modified Bcssel function of the second kind in the conven- 
tions of [22]. Divergent polynomials corresponding to subtraction constants 
in the dispersion relation have been ignored because they are annihilated by 
the Borel transform, and hence do not contribute to the Laplace sum-rules. 
The Feynman diagrams used to calculate the leading order perturbative and 
condensate contributions are illustrated in Figure 1. ^ 




Fig. 1. Feynman diagrams used to calculate the leading order perturbative and 
condensate contributions to (10). Diagram (a) corresponds to the perturbative con- 
tribution, diagram (b) is used to calculate the coefficient of (aG^), diagrams (c) 
and (d) are used to determine the contribution of {gG'^), and diagrams (e)-(g) 
are needed for the dimension-eight gluonic condensates. Additional diagrams with 
identical topologies have been omitted from the Figure. 

Although the coefficient bp, of the gluon condensate has only been calculated to 
leading order for the pseudoscalar correlation function, the one-loop coefficient 
bi is entirely determined by the leading order bo combined with the renormal- 
ization group [7] . As will be seen below, bo does not enter the Laplace sum- rules 



^ See Ref. [23] for an outline of several techniques for calculating the condensate 
contributions through the operator-product expansion and for a proof that such 
methods are equivalent. 
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Ck for A; > 0, and hence it is crucial to include hi to obtain the leading ef- 
fects from (aG^). Finally, it should be noted that the three-loop perturbative 
calculation [18] uncovered an error in the two-loop expression [19] used in the 
sum-rule analysis [7] , providing further motivation for revisiting this sum-rule 
analysis of pseudoscalar gluonium. 

The instanton contribution in (10), representing non-interacting instantons of 
size p with subsequent integration over the instanton density n{p), is expected 
to be a reasonable approximation since multi-instanton effects have been found 
to be controllable [24] . The self-dual nature of the instanton imphes that the 
instanton contributions to the pseudoscalar and scalar gluonium correlators 
only differ by an overall sign. 

The Laplace sum-rules can be constructed from (10) by using the results 
contained in [9] 



A(r) 



-2ao + ai (-6 + A^je) + a2 (tt"^ - 6 + I87E - 67 



- ^ (aG^) + Co (gG^) + dor (Os) 



- 1287r- 



2aKo{a) + {l + 2a)Ki{a) 



(16) 



Ci{t) =\ [-6ao + ai (-22 + 127^) + (Stt^ - 36 + 667^; - I87I) 



^ {aG^) - do (Os) 



256r 



'Jdnip) 



P" 



(9 - 40)0^0(0) + (3 + 7a - Aa^)Ki{a) 



(t,so) ^ j 



dt 



(17) 



Oo - 2ai log {tr) 02 (tt^ - 3 log^ (ir)) 
-h{aG^) j t'^e-'^ dt (18) 

so 

00 

+ 167r^ j dn(p)p^ j t'^+Vs {pVt) F2 (pVt) e-*" dt , 



so 



where a = p^/2r and it should be noted that the one-loop term provides the 
leading {aG"^) contribution to the Laplace sum-rules. In obtaining these ex- 
pressions, renormalization-group improvement has been implemented by set- 
ting v"^ — l/r after calculating Borel transforms [25], so that in the perturba- 
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tive corrections, a is implicitly the three-loop n/ = 3 MS running coupling 



as{y) _ 1 fh logL 1 



Pl(\og'L-\ogL-l)+P2\ (19) 



with AjYs ~ 300 MeV for three active flavours, consistent with current esti- 
mates of Q;s(Mr) [13]. 



3 Laplace Sum- Rule Analysis of Pseudoscalar Gluonium 

The instanton hquid model [26] 

dn(p)=n,5(p-pe)dp ; rie = 8 x 10"^ GeV^ , p,^-^^^L_ (21) 

will be used in the Laplace sum-rule analysis, along with the value for the 
the dimension-eight condensates obtained from vacuum saturation and the 
heavy-quark expansion [3,27] 

and the instanton estimate of the dimension-six gluon condensate [3,17] 

(gG^) ^ (0.27 GeV^) (aG^) . (23) 

The dimension-six and dimension-eight condensates are thus related to the 
gluon condensate given by the determination [28] 

aG^) = (0.07 ± 0.01) GeV^ . (24) 

In the narrow resonance(s) model 

p'^^{t) = TT 5^ F^M^S (t - W) , (25) 
resonance masses are signaled by exponential decay of the sum-rules 

^fe (t, So) = E exp (-MV) . (26) 

A bound on the mass M of the hghtest state can then be obtained from the 
ratio 

¥f\= lim |4^>M^ . (27) 
Lo{t) 5o->oo (r. So) 
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This bound is quite robust since it does not depend on the QCD continuum 
approximation and does not require dominance from the hghtest state (see 
e.g. [9]). Figure 2 displays the ratio (27), resulting in a bound on the lightest 
pseudoscalar gluonium state of M < 3.1 GeV. 




Fig. 2. The ratio ^ Lxj as a function of r for central values of the QCD param- 
eters. 

As illustrated in Figures 3 and 4, the rise past the minimum in Figure 2 is 
actually a manifestation of a zero which occurs in Lq{t) prior to the first 
zero of vCi(t), resulting in a singularity of the ratio (27). Negative values of 
£fe(r) are inconsistent with (5) and positivity of the spectral function pit). 
This unphysical behaviour of the sum-rules can be traced to the instanton 
contributions in the low-energy (large r) regime. 




Fig. 3. The sum-rule £o (t") = liniso^oo -^^o (''") ^o) as a function of r. Central values 

of the QCD parameters have been used. 

The unphysical low-energy behaviour of the sum-rules does not necessarily 
invalidate the bounds obtained from Figure 2 because the minimum occurs in 
the perturbative (high-energy, low r) regime, but it does present difficulties 
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Fig. 4. The sum-rule £i (r) = lim^jg^oo Ri (t, sq) as a function of r. Central values 
of the QCD parameters have been used. 

in the extraction of an optimized mass estimate. Within the single narrow 
resonance model, (26) results in 



Ro (t, So) 



= . (28) 



In a typical sum-rule analysis [17], one decreases sq until a stable {i.e. nearly 
r-independent) ratio (28) is obtained. Figure 5 shows that the ratio does sta- 
bilize at approximately 2.7 GeV, but the small r values associated with this 
stability region are problematic because dominance of the lightest resonance 
cannot be assured and uncertainties from the continuum approximation be- 
come significant. However, because (27) has general validity beyond these 
approximations, the M < 3.1 GeV bound remains valid. 



ID 




□ D D1 D2 as D4 DE 06 



Fig. 5. Sum-rule ratios ^ Ri (r, sq) /Rq (r, sq) ranging from sq = oo {i.e. yZi/Zo) 
in the top curve to sq = 8 GeV^ in the bottom curve. Central values of the QCD 

parameters have been used. 

Thus instanton effects prevent the extraction of a reliable optimized Laplace 
sum-rule mass prediction for pseudoscalar gluonium, a difficulty that has also 
been observed in the configuration-space instanton analysis [24] . It is interest- 
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ing that the resonance signal in each case is exponential decay associated with 
the resonance mass combined with r in the Laplace sum-rules (26) or with 
Euclideanized time in the configuration-space approach [24] . The next Section 
will review the formulation of Gaussian sum-rules which will be seen to pro- 
vide a fundamentally-different weighting of the hadronic spectral function and 
QCD contributions, obviating the problems associated with an optimization 
of the Laplace sum-rule analysis. 



4 Gaussian Sum-Rules for Pseudoscaleir Gluonium 

The simplest Gaussian sum-rule (GSR) [29] 




relates the QCD prediction on the left-hand side of (29) to the hadronic spec- 
tral function p{t) smeared over the energy range s — 2y^ ^ i s + 2y/T, 
representing an energy interval for quark-hadron duality. ^ An interesting as- 
pect of the GSR is that the duality interval is actually constrained by QCD. 
A lower bound on this duality scale r necessarily exists because the QCD pre- 
diction has renormalization-group properties that reference running quantities 
the the energy scale — ^JT [29,30]. Thus it is not possible to achieve the 
formal r ^ limit where complete knowledge of the spectral function could 
be obtained via 

lim G (s, t) = -p (s) , s>tQ . (30) 

However, there is no theoretical constraint on the quantity s representing the 
peak of the Gaussian kernel appearing in (29). Thus the s dependence of the 
QCD prediction G (s, r) probes the behaviour of the smeared spectral function, 
reproducing the essential features of the spectral function. In particular, as s 
passes through t values corresponding to resonance peaks, the Gaussian kernel 
in (29) reaches its maximum value, implying that Gaussian sum-rules weight 
excited and ground states equally. 

The QCD prediction for the GSR is obtained from the correlation function 
through 

;iB{ "(-^-^)r/(-^+-^) } , (31) 

^ The quantity r in the GSR should not be confused with the similar quantity 
appearing in the Laplace sum-rules. In particular, the mass dimension is different 
in each usage. 
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where the Borel transform B has been redefined as 



B = Um 

A^/N=4t 



-A 



2\N 



d 



JV 



(32) 



The GSR (29) then results from combining (31) with the dispersion relation 
(2) [29]. 

A connection between Gaussian and finite-energy sum-rules can be established 
through the diffusion equation satisfied by the GSR 



d^G{s,T) dG{s,T) 



dr 



(33) 



In particular, the rcsonancc(s) plus continuum model (6), when p^"''^{t) is 
evolved through the diffusion equation (33), only reproduces the QCD predic- 
tion at large energies (r large) if the resonance and continuum contributions 
are balanced through the finite-energy sum-rule [29] 



1 ^° 

F{so) = - f p^'"'(t)&t 

TT J 

to 



(34) 



Within the resonance (s) plus continuum model (6), the continuum contribu- 
tion to the GSR is determined by QCD 



G'^°"*(s,r,so) = 



1 



TTT Jso 



exp 



At 



TT 



-lmn'2^^(i)di , 



(35) 



and is thus combined with (<§, t) to give the total QCD contribution 

G^'" {s, T, so) = G {s, t) - G^' {s, T, so) , (36) 
resulting in the final relation between the QCD and hadronic sides of the GSR. 



G^^''{s,T,So) 



Jto 



exp 



t)' 



At 



TT 



(37) 



Comparison of (37) and (8) reveals that the GSR provides a fundamentally 
different weighting of the hadronic spectral function than the Laplace sum- 
rule. In particular, the Laplace sum-rules always emphasize the low-energy 
region, while this aspect is controlled in the GSR by the parameter s. 



The expression for the QCD prediction for the pseudoscalar gluonium GSR 
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can then be constructed using the results of [16]. 



1 r^o 
V47rT Jo 



At 



(ao - 7r^a2) + 2ai log ( — 



+3«2log'(^ 



1 - - 



— --hi{aG ) / exp 
s/A-KT Jo 



4r 



dt 




Co (Oe) - ^ (Os) 



TTT 



dn(p)p' / 
Jo 



-js-tf 
At 



(38) 



As mentioned earlier, renormalization-group improvement necessitates the re- 
placement z/^ = ^/T in (38) [29,30]. As with the Laplace sum-rules, the one-loop 
term provides the leading (aG"^) contribution to the GSR. For the dimension- 
six and -eight non-perturbative QCD condensate contributions in (38), it is 
easily seen that the non-perturbative corrections are exponentially suppressed 
for large s. Since s represents the location of the Gaussian peak on the phe- 
nomenological side of the sum-rule, the non-pcrturbative corrections are most 
important in the low-energy region, as anticipated by the role of QCD conden- 
sates in relation to the vacuum properties of QCD. This explicit low-energy 
role of the QCD condensates clearly exhibited for the Gaussian sum-rules is 
obscured in the Laplace sum-rules. 

Integrating both sides of (37) reveals that the overall normalization of the 
above equation is related to the finite-energy sum-rule [30] 



00 1 °° 

J G'i'^is, T, So) = - / dt 

-00 to 



(39) 



Thus the diffusion equation analysis [16] relates the normalization of the GSR 
to the finite-energy sum-rules. Information independent of this relation is thus 
extracted from the normalized GSR [30] 



iV'"='(5,r,so) 



Gi^''{s,T, So) 



(40) 



mo (r. So) 

00 

,,(r,So)= / W^"(s,r,So) ds , n = 0,l,2,... , (41) 
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which is related to the hadronic spectral function via 



/exp 



4t 



p'"^{t) dt 



to 



(42) 



5 Gaussian Sum-Rule Analysis of Pseudoscalar Gluonium 



The techniques for analyzing the GSRs were initially developed in [16,30]. In 
the single narrow resonance model, the normalized GSR (42) becomes 



7V^'='^(s,T,So) = ^=exp 



MY 



At 



(43) 



Deviations from the narrow- width limit arc proportional to M^F^/r, so this 
narrow-width model may be a good numerical approximation. Phenomenolog- 
ical analysis of the single narrow resonance model proceeds from the observa- 
tion that, as a function of s, the phenomenological side of (43) has a maximum 
value (peak) at s = independent of the value of r. The value of sq is then 
optimized by minimizing the r dependence of the s peak position of the QCD 
prediction, and the resulting r-averaged s peak position leads to a prediction 
of the resonance mass 



For the central values of the QCD parameters, minimizing the s peak motion in 
the region 2 GeV^ < r < 4 GcV^ results in M = 2.70 GeV for ^ = 3.3 GeV, 
which is remarkably similar to the mass scale resulting from the Laplace sum- 
rule stability analysis in Figure 5.^ Although an optimized Laplace sum- 
rule analysis was not possible, a consistent scenario is emerging from the two 
approaches. 



Figure 6 compares the theoretical and phenomenological sides of the normal- 
ized GSR (42). It is evident from this comparison that the single narrow res- 
onance model is an inadequate description of the spectral function predicted 
by QCD. Figure 6 also reveals a region of small s where the theoretical GSR 
is negative, inconsistent with a positive spectral function through (42). This 
effect can again be traced to the instanton contributions, and since s corre- 
sponds to the peak of the Gaussian kernel, it is clear that this is a low-energy 
non-perturbative effect that is safely isolated from the peak of the theoretical 
contribution which enters the optimization analysis. 



The range of r is chosen to have acceptable convergence of the perturbative series 
while maintaining a resolution consistent with typical hadronic scales. 
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By construction, all the curves in Figure 6 are normalized to unit area, so 
the QCD contributions which underestimate the peak are necessarily broader 
than those of the single narrow resonance model. In particular, moment 
combinations (41) associated with (43) result in 



mo \moJ 



(44) 



Thus the narrow-resonance width of 2t is insufficient to provide agreement 
with QCD, indicative of distributed resonance strength that can be resolved 
by the GSR. 




s (GeV ) 



Fig. 6. Comparison of the QCD and single narrow resonance model contributions to 
the normalized GSRs as a function of s for the optimized parameters M = 2.70 GeV 
and ^/so = 3.3 GeV for r = 2 GeV^ (upper pair of curves) and for r = 4 GeV^ (lower 
pair of curves). In each case, the peak of the theoretical (QCD) prediction lies below 
that of the single-resonance model. 



A simple toy model that illustrates the effect of resonance widths in the 
broadening of the phenomenological side of the normalized GSR is a unit- 
area "square pulse" which could describe a broad structureless feature of the 
spectral function 



e {t-M^ + my) -e{t-M^ -my 



(45) 



which leads to the following normahzed Gaussian sum-rule [16] 



AMY 



erf 



's-M^ + MY^ 



— erf 



's-M^- MY\ 



(46) 



^ This point is explicitly evident in the tails where QCD overestimates the single 
narrow resonance model. 
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The second-order moment combinations resulting from this sum-rule are [16] 



mo Vmo/ 3 



(47) 



and hence resonance widths broaden the phenomenological side of the nor- 
malized GSR as would be expected intuitively. 

For detailed analysis, the following Gaussian resonance and a skewed-Gaussian 
resonance models will be employed since they are numerically simpler to an- 
alyze than a Breit-Wigner resonance shape. 



2\2 



p{t) ~ exp 
p{t) ~ exp 



{t - M ^) 

{t-M 



2^2 



(48) 
(49) 



The quantity 7 can be related to a Breit-Wigner width F by equating the half- 
widths, resulting in F = y/2\og2^/M. For the Gaussian model, the resulting 
normalized GSR is [16] 



l + erff^2^±^!^ 



exp 



'2(72 + 2t) 



(50) 



For the skewed Gaussian model, the resulting normalized GSR is 



7V^^'^(S,T, So) 



= exp 



2n2 1 



"2(72 + 2t) 



1 + erf 



/ S72 + 2MV y 



X s^7^ + 4s7Vm2 + 4t2M^ + 2t7^ + 4x^7^ 



(51) 



+ 



27 (s7^ + 2rM2) (72 + 2t) exp 



s^7^ + 2rM4 
4x7^ 



where 



20F(72+2r)5/2 



M27^ exp (-^) + (M^ + 7=^) [1 + erf {^) 



(52) 



The inclusion of resonance widths implies that the s peak position of the phe- 
nomenological side of the sum-rules can develop r dependence, complicating 

the optimization of Sq. The direct approach of determining all parameters from 
the best fit between the s dependence of the two sides of (50) and (51) over 
the range 2 GeV < r < 4 GeV^ is facilitated by an initial estimate of Sq. An 
effective estimate of the optimized sq is obtained from the (approximate) r 
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dependence of the s peak which has the general behaviour [16,30] 

be 

Speak (r, So) = a H \-— . (53) 

T 

Analysis of how the s peak "drifts" with r in comparison with the behaviour 
(53) provides an estimate of Sq which is found to be surprisingly close to the 
true fitted value, facilitating numerical analysis of the multi-dimensional fit. 

For the central values of the QCD parameters, the resulting optimized pa- 
rameters are M — 2.7 GeV, F = 0.50 GeV in the Gaussian resonance model, 
and M = 2.6 GeV, F = 0.53 GeV in the skewed Gaussian resonance model, 
with = 3.2 GeV in each case. The agreement between the QCD and phe- 
nomenological sides of the normalized GSR for the skewed Gaussian model 
is shown in Figure 7. Comparison of Figures 6 and 7 (which have identical 
scales) reveals that inclusion of a resonance width provides a significant im- 
provement in the agreement between the QCD and phenomenological sides of 
the sum-rule without a significant change in the predicted resonance mass. 



□ ZD 




Fig. 7. Comparison of the QCD and skewed-Gaussian resonance model contri- 
butions to the normalized GSR as a function of s for the optimized parameters 
M = 2.6 GeV, T = 0.53 GeV, and ^ = 3.2 GeV for r = 2 GcV (upper pair of 
curves) and for r = 4GeV^ (lower pair of curves). For the upper pair of curves 
the peak of the theoretical (QCD) prediction lies above that of the single-resonance 
model, but the situation is reversed for the lower pair of curves. 

The un-skewed Gaussian resonance model results in an agreement between 
the QCD and phenomenological sides of the normalized GSR which is quali- 
tatively and quantitatively indistinguishable from Figure 7. Unlike the single 
narrow resonance model where theory consistently underestimated the phe- 
nomenological peak, there is no evident pattern in the disagreement between 
QCD and the skewed-Gaussian model. It is possible that the slight disagree- 
ment that remains between QCD and the phenomenological model can be 
explained by the low-energy region where G'^'^'^ is negative, since this reduces 
the quantity mo for r = 2 GeV*^, and hence raises the value of A^^^'^. Thus 
there is no clear motivation for considering more complex resonance models. 
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The effect of varying tlie QCD instanton parameters by 15% and the gluon 
condensate within the range (24) in both the Gaussian and skewed- Gaussian 
models provides a final determination of the mass and width of pseudoscalar 
gluonium 

M = (2.65 ± 0.33) GeV , r< 540MeV . (54) 
An interesting aspect of this analysis is that smaller instanton sizes are essen- 
tially consistent with a narrow resonance model. 



6 Conclusions 



Instanton and higher-loop perturbative effects have been incorporated into 
the Laplace and Gaussian sum-rules for pseudoscalar gluonium. Although the 
Gaussian sum-rules are found to be superior to the Laplace sum-rules in ob- 
taining an optimized prediction of the pseudoscalar gluonium properties, the 
resulting mass scales in each analysis are remarkably consistent, providing cor- 
roborating evidence in support of our final mass prediction of M ^ 2.7 GeV. 
Our mass estimate is consistent with the benchmark lattice scalar glueball 
mass of approximately 1.6 GeV and a pseudoscalar to scalar mass ratio of 
approximately 1.5 [1]. 

The analysis presented in this paper has not pursued the interesting possibility 
of mixing between gluonium and quark mesons [31], although in principle any 
hadronic state that overlaps with the pseudoscalar gluonium interpolating field 
would be probed by the correlator (1) which underlies our entire analysis. 

This possibility of other resonances coupled to the pseudoscalar gluonic cur- 
rent was explored through an extension of the Section 5 analysis to two narrow 
resonances tends to lead to nearly-degenerate states near 2.7 GeV, suggesting 
that the 2.7 GeV state is more strongly coupled to the pseudoscalar gluonic 
operator than other pseudoscalar states such as the Vj . This result is not unex- 
pected, since the coupling of the 2.7 GeV state |P) to the pseudoscalar gluonic 
current extracted from the peak value of G^'^'^ (s, r, sq) with the optimized value 
of So is 



aGG 



P) 



20GeV3 



(55) 



By comparison, the couplings of 77, 77' to the pseudoscalar gluonic operators 
are [32] 



aGG 



v') 



0.9 GeV^ 



{O 



aGG 



v) 



0.1 GeV^ 



(56) 



and hence the relative strength of the \ri') compared with \P) obtained from 
the ratios of the squares of (55) and (56) is approximately 0.2%, and the rela- 
tive strength of the \ri) is approximately 0.02%. This result is upheld by a fit of 
the relative strength of the rj' and a Gaussian resonance to the normalized GSR 
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results in a relative strength of less than 1%. Thus the 2.7 GeV state result- 
ing from our analysis appears to be dominantly-coupled to the pseudoscalar 
gluonic operators. 
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